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Numerical Solution of the Frequency Equations for the 
Flexural Vibration of Cylindrical Rods 

W. E. Tefft 

(June 28 7 I960) 

A numerical solution of Pickett's frequency equations for the flexural vibration of 
cylindrical rods, based on the three-dimensional differential equations of elasticity, 1ms 
been obtained on the IBM 704 computer. The results are presented in the form of tables 
of correction factors to be applied to the thin rod approximation for the fundamental 
flexural resonance frequency and the first two overtones. These results provide an accurate 
means of calculating Young's modulus from the density, dimensions, and resonance fre- 
quencies of cylindrical rods having diameter-to-length ratios as high as 0.(>. 



1. Introduction 

Improved experimental techniques recently have 
made it possible to determine the elastic properties 
of materials to a high degree of precision by dynamic 
resonance methods. Factors which have limited 
accurate measurements are: Convenient experimental 
techniques are not well developed theoretically, and 
conversely, techniques based on well analyzed 
theories are experimentally difficult. A recent ex- 
perimental investigation [1] ' showed that Pickett's 
[2] approximate solution of the flexure problem, 
based on Goen's [3] solution of Timoshenko's [4] 
equations may be in error by a feu percent. In a 
later paper, Pickett [5] gave an accurate set of 
frequency equations based on the three-dimensional 
differential equations of elasticity. The purpose of 
the present paper is to present the results of a 
numerical solution of these more accurate 4 equations. 
These numerical solutions have put the accurate 
theory of flexural vibration of cylinders into a more 
usable form. Elimination of theoretical errors, for 
this particular problem, has thus been achieved. 

2. Background 

According to the elementary, one-dimensional 
theory, Young's modulus, E, may be computed from 
the flexural resonance frequency, /, of a cylinder of 
density p, length I, and diameter d, from the equa- 
tion [3] 
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E= 



mid 2 
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where f n is the resonance frequency of the n th mode 
of vibration and the m n are solutions of the equation 

cos m cosh m = 1 . 

for the first three modes of vibration 



The values of m 
are [6] 



and 



m! = 4.730041, 
m 2 = 7.853205, 
m 3 = 10.995608. 



(3) 
(4) 

(5) 



1 Figures in brackets indicate the literature references at the end of this paper 



Timoshenko [4] added shear and rotatory inertia 
corrections to the differential equation, and (mens [3] 
reduced the resulting frequency equation to the form 
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where T n is a function of d/l and Poisson's ratio, 
ijl// 1 ,, approaches one as d/l approaches zero, for 
all values of // and p. 

Pickett [2] made 4 numerical computations of the 4 
T n , and gave approximate equations for this com- 
putation. These values are not as accurate as is 
desirable, in view of the experimental precision now 
attainable in dynamic resonance measurements. 
The basis for the more accurate values presented 
in this paper was given by Pickett [5] in his solution 
based upon the three-dimensional differential equa- 
tions of elasticity. 

A brief verbal sketch of Pickett's derivation may 
be helpful in understanding the meaning of some of 
the following equations. He begins by rioting that 
there are three fundamental solutions for elastic 
waves in a cylinder. Each is characterized by the 
angular frequency p and the wavelength 2ira/y, 
where a is the radius of the cylinder. These three 
solutions can be combined so that their sum satisfies 
the boundary conditions on the curved surface of the 
cylinder. This gives a relationship between p and 
7 expressed by Pickett's [5] equation 12. For each 
value of p there are two values of 7, denoted by 7 
and 7', which satisfy this equation. One of these, 
y', corresponds to a sum which represents a traveling 
wave. The other, 7, corresponds to a sum in which 
the "wavelength" is an imaginary number, i.e., the 
displacement is an exponential, rather than a sinu- 
soidal function of distance along the rod. Pickett 
combines both of these sums to get a solution for 
displacements which approximately satisfies the 
boundary conditions on the flat ends of the cylinder. 
These boundary conditions lead to his equation 21, 
the frequency condition. The present paper is 
concerned with the simultaneous numerical solution 
of these equations. 
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3. Theory and Calculations 

Pickett's [5] results were reduced to the form of 
equations (7) through (12) below, by a method similar 
to that used by Hudson [7]. Pickett's equation 12 is 
an expression of the requirement that the three radial 
components of stress be zero at the surface of the 
cylinder. This equation was reduced to the form of 
eq (7) by the aid of definition (10) and a suitable 
manipulation of the resulting determinant to reduce 
it from third order to second order. 

2K^(K) + (K 2 -y 2 )[2iP(K) +K 2 \ -2y 2 [^(K) +K 2 

2K 2 +(h) + (K 2 +^ 2 ) (K 2 -y 2 )*(K) + {K 2 +y 2 , 

= 0, (7) 
where h and K are defined by 
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where 

jit=Poisson's ratio, 
E= Young's modulus, 
p= density, 

a=d/2, d being diameter of cylinder, 
P=2tJ= angular frequency of flexural vibration, 
7= parameter determined from the solution of 

these equations, 
[//=& combination of Bessel functions of the first 

kind, orders zero and one, and is defined by 
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where the upper terms in braces refer to the odd 
numbered and the lower terms to the even numbered 
modes of vibration, then equation 21 of [5] becomes 



F(y)=Ftf). 



(12) 



The problem of determining the flexural resonance 
frequencies is now solved in principle. For a given 
length, diameter, Young's modulus, and Poisson's 
ratio, one must solve (7) and (12) simultaneously for 
p (and incidentally for 7 and y f ). However, these 
are transcendental equations and involve the fre- 
quency in such a complex way that numerical meth- 
ods must be used. 

A program for obtaining a numerical solution of 
the simultaneous set of transcendental eqs (7) and 
(12), using (11), on the IBM 704 computer was 
written by the Bureau's computer section. The 
numerical solution was made in the manner sug- 
gested by Pickett [5]. Appropriate values of pa 2 p 2 /E 
and /x were chosen. The two real values of y 2 satisfy- 
ing eq (7) were determined by successive approxi- 
mations. These values of y 2 were then substituted 
in eq (11) and the values of l/2a satisfying eq (12) 
for each of the first three modes of vibration were 
determined by successive approximation. The cor- 
rection factors, T n , were then computed from the 
relation 



a p q2 P 2 
E 



(13) 



Twenty-nine values of pa 2 p 2 /E were chosen, rang- 
ing in value from LOXIO" 6 to 0.3958. For each of 
these values of pa 2 p 2 /E, eleven values of p were 
chosen, ranging from zero to 0.50 in steps of 0.05. 
There were thus 319 computed points for each of 
the first three modes of vibration. 

The values of T lf T 2 , and T z were then computed 
for equal intervals of 2a /I by a sixth-order Aitken's [8] 
interpolation method. The overall error in these 
computations was less than three parts in 10 6 for 
values of pa 2 p 2 /E greater than about 0.0007, corre- 
sponding to d/l — 0.1 for the fundamental mode. For 
lower values of pa 2 p 2 /E f however, machine round-off 
error became significant. In order to obtain accu- 
rate values of T for pa 2 p 2 /E<^0. 0007 , which corre- 
sponds to long thin rods, it was necessary to make 
a first-order approximation to the set of eqs (7), (11), 
and (12). This was accomplished as follows: 

If we let 

(l+M)(l-2„) 2(1+M ) = jB ,and^=P, 

1 — p Hi 

then, from eqs. (8) and (9), 

h 2 =AP-y 2 , K 2 =BP-y 2 , (14) 

and from eq (10) we obtain 



m- 



AP—fT, . AP-y 2 . (AP-y 2 Y 2 
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and 



BP-y 2 F, , BP-y 2 . (BP-y 2 ) 2 



I 



24 ' 384 
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If eqs (14), (15), and (16) are substituted in eq (7), 
terms of order y 2 cancel (terms of the form y 2 p n are 
higher order) . Inspection of the lowest order terms 
present shows that 



7 2 ~±2VP. 



(17) 



But obtaining a solution to (12) requires that the 
magnitude of y 2 and y' 2 be different. We must, 
therefore, make a second-order approximation to 
eq (7). Also, it is obvious from eq (17) that y 2 is of 
order -yjP, not of order P; therefore, we must keep 
the first three terms of (15) and (16) for this approxi- 
mation. When this is done, and algebraic simplifica- 
tions are made, we arrive at 



7 4 -4P-^(5 + 2m)Pt 2 



12(1+m) 



=0, (18) 



where the last two terms are of higher order than the 
first two. In order to solve (18) for the two real 
values of y 2 , assume 



T*=Cm ! P+DP, 



(19) 



where C and D are functions of p. to be determined. 
The substitution of (19) into (18) results in a poly- 
nomial in ^P. If we neglect terms of order higher 
than p d/2 , and equate the coefficients of each of 
the remaining powers of P separately to zero, we 
arrive at the following values for C and D. 



and 

Therefore, 

7 
and 



C=±2 

10+15/x + V 



Z>=- 



6(1 +m) 



=+ 2 Vt + v 6(i+ M ) jve-) m 

7 - ~ 2 V^rH 6(i+m) ; \te-/ (21) 

Equations (20) and (21) were then substituted in 
(11) and (12), and the functions cot 7^/2a, cot y'l/2a, 
tan yl/2a, and tan y'l/2a were approximated in the 
following manner: since 7 2 <C0, y is imaginary, so 
the following relations hold identically: 



and 



cot ~-=— i coth ~~ 
2a 2a 



tan 7r=i tanh ^— 
2a 2a 



Furthermore, it is known from the elementary solu- 
tion that 



lim cot 7^-=lim ( —coth -Q - )=cot 
2a 2a 2a \ 2a J 



for odd-numbered modes, and 

lim tan +- =lim tanh -Q— =tan — 
2« ^ 2a 2a , 2a 2 

for even-numbered modes. We therefore make a 
first-order Taylor's series expansion of these functions 
about the point m/2, giving 

y'l m 



cot 



2a ' 



-cot 



m 2a 



sur 



m 



and similar approximations for the other three 
functions. 

The quantity pay 2 / 1 ^ was eliminated by means of 
(13), and the various powers of T which appeared 
were eliminated by means of a Taylor's series expan- 
sion about the point d/l=0, T=l. This expansion 
can be shown to contain no terms of order d/l, so the 
first correction term in the expansion is proportional 
to (d/l) 2 . 

With these substitutions and further algebraic and 
trigonometric manipulation, we arrive at the follow- 
ing equation 



T«l + 



[( 



m 2 (sin m+sinh m)\ /'l()+15/x+ 4m 
sinh m— sin : 

/m sin m sinh m 
\ sin m— sinh m 



48(1+m) 
'2-3/*- V 



) 



n\ / 2-3m-V \1 /A 2 
<i)\ 24(l + /x) )\\l) 



where the upper (plus) sign refers to odd-numbered 
modes and the lower (minus) sign refers to even- 
numbered modes of vibration. 

Inserting the values of m from eq (3), (4), and (5) 
we obtain 



2\« l+4.88669( 



Vl + 1.26225/z+0.20980/x 2 \ (d\ 2 

\~ ~r+jr -)\t) 

T 2 ~i+n.,2 ^+ 1 ^ 



r 3 « 1 + 26.10281 



(■ 



1 + 1 .39469^+0.31575^ 



)©' 



(22) 



(23) 



(24) 



A plot of T versus (d/l) 2 is nearly a straight line. 
This indicates that the first term (of order (d/l) 2 ) in 
the Taylor's series expansion of T gives most of the 
contribution over the range of interest. Thus, the 
function 

varies only very slowly with (d/l) 2 



5627-5.9 — 60h- 
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The method for obtaining accurate values of T in 
the region (d/Z) 2 <<l is now apparent. A plot of 
H versus (d/l) 2 is made on a highly expanded scale, 
for each of the first three modes of vibration and for 
each value of jjl, with the H intercept being obtained 
in each case from the appropriate one of eq (22), 
(23), or (24) and the other points being taken from 
the computer solution. When plotted in this manner, 
the first few points from the computer solution did 
not form a smooth curve. This scattering of points 
was an indication of significant round-off error at 
these low values of (d/l) 2 . These points were there- 
fore ignored, and the smooth curve, obtained at 
higher values of (d/l) 2 , was extended to the chosen 
intercept. Finally, the values of H for equally 
spaced intervals of (d/l) were read from this smoothed 
curve, and the values of T were computed from eq 
(25). 

4. Results and Discussion 

Tables 1,2, and 3 give T as a function of d/l and 
fi, for the first, second, and third modes respectively. 
These tables are accurate to about three parts in 10 6 . 
In order to calculate Young's modulus from measure- 
ments of the flexural resonance frequencies of 
cylinders, the following equations, based on eq (6), 
are used: 

£'=1.261886^^ Tit 



E-- 



= 0.1660703^? T 2 t 



and 



#=0.04321184 



d 2 



t8* 



for the first, second, and third modes respectively. 

A linear interpolation of T versus (d/l) 2 or T 
versus p can be made with an accuracy of better than 
one part in 10 4 over the range of the tables. To 
obtain higher accuracy, Aitken's interpolation 
method should be used. 

In order to obtain maximum accuracy from this 
method of Young's modulus determination, it is 
recommended that, whenever feasible, resonance 
frequencies of the fundamental and first two over- 
tones be measured, and the values of E be computed 
for each. An average of these values of E should be 
more accurate than any of them alone. The reason 
for this is that small irregularities in the diameter of 
a specimen, which can never be eliminated com- 
pletely, will affect each mode differently; i.e., some 
modes will have their frequencies shifted upward, 
and some downward. The magnitude of the dis- 
agreement among these values of E indicates the 
overall uniformity of the dimensions of the specimen. 

The correctness of these computations was checked 
experimentally, by a method which will be described 
in a subsequent publication. Table 4 gives a com- 
parison of experimental and theoretical values of T 
for two of the most accurately machined specimens. 
The discrepancy between theory and experiment is 
seen to be of the order of 0.03 percent, except for 
T 3 . The experimental value of T 3 was difficult to 



Table 1. The correction factor for the fundamental mode of flexwal vibration of cylinders, T\, as a function of the diameter to 

length ratio, d/l, and Poisson's ratio, /*. 

The values of T\ for d/Z<0.08 were obtained by a graphical interpolation as explained in the text. All other values of T\ were obtained from the computer solution. 



d/T\ m 


0.00 


0.05 


0.10 


0.15 


0.20 


0.25 


0.30 


0.35 


0.40 


0.45 


0.50 


0.00 


1.000000 


1.000000 


1.000000 


1. 000000 


1. 000000 


1.000000 


1. 000000 


1. 000000 


1. 000000 


1. 000000 


1. 000000 


.02 


1.001954 


1.001979 


1.002001 


1. 002029 


1.002053 


1.002077 


1. 002100 


1.002124 


1. 002147 


1. 002170 


1. 002193 


.04 


1.007804 


1. 007905 


1.008004 


1. 008102 


1. 008199 


1. 008295 


1.008388 


1. 008482 


1. 008575 


1.008666 


1. 008757 


.06 


1. 017522 


1. 017748 


1.017968 


1. 018186 


1. 018405 


1. 018619 


1. 018826 


1. 019038 


1. 019245 


1. 019450 


1. 019653 


.08 


1. 031064 


1. 031461 


1.031848 


1. 032233 


1.032618 


1.032994 


1.033360 


1. 033733 


1. 034096 


1. 034459 


1. 034818 


.10 


1. 048367 


1. 048983 


1. 049580 


1. 050174 


1. 050765 


1. 051344 


1.051916 


1. 052484 


1. 053050 


1.053610 


1. 054170 


.12 


1.069357 


1. 070225 


1.071077 


1. 071920 


1. 072753 


1.073577 


1. 074393 


1. 075202 


1. 076008 


1. 076808 


1.077604 


.14 


1.093950 


1.095111 


1. 096256 


1. 097378 


1. 098495 


1. 099599 


1. 100694 


1.101782 


1. 102864 


1. 103941 


1. 105015 


.16 


1. 122061 


1. 123544 


1. 125007 


1. 126452 


1. 127884 


1. 129302 


1. 130711 


1. 132113 


1. 133509 


1. 134901 


1. 136288 


.18 


1. 153592 


1. 155430 


1. 157245 


1. 159039 


1. 160817 


1. 162582 


1. 164337 


1. 166086 


1. 167830 


1. 169569 


1. 171308 


.20 


1. 188458 


1. 190676 


1. 192869 


1. 195038 


1. 197191 


1.199332 


1. 201464 


1. 203591 


1. 205714 


1.207838 


1. 209961 


.22 


1. 226567 


1.229189 


1. 231781 


1. 234351 


1. 236906 


1.239449 


1.241986 


1. 244521 


1. 247057 


1. 249596 


1. 252141 


.24 


1. 267838 


1. 270881 


1. 273894 


1. 276886 


1. 279865 


1.282836 


1. 285807 


1. 288779 


1. 291757 


1.294746 


1. 297747 


.26 


1.312189 


1. 315668 


1.319120 


1. 322555 


1. 325980 


1. 329403 


1. 332832 


1. 336270 


1. 339723 


1. 343194 


1. 346687 


.28 


1.359546 


1. 363477 


1. 367381 


1. 371276 


1. 375167 


1. 379065 


1. 382977 


1. 386911 


1. 390869 


1. 394857 


1. 398879 


.30 


1. 409842 


1. 414233 


1. 418606 


1. 422974 


1. 427352 


1. 431747 


1. 436169 


1. 440625 


1. 445121 


1. 449661 


1. 454250 


.32 


1. 463014 


1. 467873 


1. 472724 


1. 477584 


1. 482465 


1. 487380 


1.492337 


1.497345 


1. 502411 


1. 507540 


1. 512737 


.34 


1. 519006 


1.524340 


1. 529680 


1. 535043 


1. 540446 


1. 545901 


1. 551420 


1. 557012 


1. 562682 


1. 568439 


1. 574286 


.36 


1. 577766 


1.583580 


1. 589417 


1. 595298 


1.601240 


1. 607259 


1.613367 


1. 619572 


1. 625884 


1. 632309 


1. 638853 


.38 


1. 639253 


1. 645548 


1. 651888 


1. 658300 


1. 664800 


1. 671405 


1. 678129 


1. 684983 


1. 691975 


1.699113 


1. 706401 


.40 


1. 703422 


1.710202 


1. 717054 


1. 724007 


1. 731082 


1. 738298 


1. 745669 


1. 753207 


1. 760921 


1. 768818 


1. 776904 


.42 


1. 770243 


1. 777505 


1. 784874 


1. 792382 


1. 800052 


1. 807904 


1. 815954 


1. 824214 


1. 832694 


1.841401 


1. 850342 


.44 


1. 839681 


1. 847424 


1. 855317 


1. 863393 


1.871677 


1. 880193 


1. 888956 


1. 897979 


1. 907273 


1.916846 


1.926702 


.46 


1.911712 


1. 919934 


1. 928355 


1. 937012 


1. 945932 


1. 955140 


1. 964653 


1. 974485 


1.984646 


1. 995143 


2. 005980 


.48 


1.986311 


1. 995009 


2. 003964 


2. 013216 


2. 022795 


2. 032727 


2. 043030 


2. 053718 


2. 064802 


2. 076288 


2. 088178 


.50 


2. 063458 


2. 072628 


2. 082123 


2. 091985 


2. 102247 


2.112937 


2. 124075 


2. 135673 


2. 147741 


2. 160285 


2. 173305 


.52 


2. 143136 


2. 152775 


2. 162814 


2. 173303 


2. 184276 


2. 195762 


2. 207781 


2. 220345 


2. 233465 


2. 247142 


2. 261377 


.54 


2. 225331 


2. 235431 


2. 246024 


2. 257157 


2. 268871 


2. 281194 


2. 294146 


2. 307741 


2. 321984 


2. 336876 


2. 352415 


.56 


2. 310030 


2. 320588 


2. 331739 


2. 343539 


2. 356026 


2. 369231 


2. 383174 


2. 397866 


2. 413311 


2. 429508 


2. 446448 


.58 


2. 397222 


2. 408232 


2. 419951 


2. 432439 


2. 445736 


2. 459873 


2. 474869 


2. 490733 


2. 507467 


2. 525064 


2. 543511 


.60 


2. 486902 


2. 498357 


2. 510654 


2. 523855 


2. 538002 


2. 553126 


2. 569244 


2. 586362 


2. 604478 


2. 623578 


2. 643647 
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Table 2. The correction facto? for the first overtone of flexural vibration of cylinders, T 2 , as a function of the diameter to length 

ratio, d/l, and Po.sson's ratio, n. 

The values of Ti for rf/Z<0.05 were obtained by a graphical interpolation, as explained in the text. All other values were obtained from the computer solution. 



d/T\ m 


0.00 


0.05 


0.10 


0.15 


0.20 


0.25 


0.30 


0.35 


0.40 


0.45 


0.50 


0.00 


1.000000 


1.000000 


1.000000 


1.000000 


1 . 000000 


1 . 000000 


1.000000 


1.000000 


1.000000 


1.000000 


1.000000 


.01 


1.001352 


1.001376 


1.001399 


1.001422 


1.001445 


1.001467 


1.001489 


1.001511 


1.001533 


1.001554 


1.001576 


.02 


1.005404 


1.005499 


1.005592 


1.005683 


1 . 005774 


1.005863 


1.005952 


1 . 006039 


L. 006126 


1.006212 


1 . 006297 


.03 


1.012144 


1.012355 


1.012564 


1.012770 


1.012973 


1.013174 


1.013373 


1.013569 


1.013763 


1.013057 


1.014149 


.04 


1.021552 


1.021926 


1.022295 


1.022660 


1.023019 


1.023376 


1.023728 


1.024077 


1.024421 


1.024763 


1.025103 


.05 


1.033605 


1.034188 


1.034760 


1.035326 


1.035885 


1.036440 


1.036987 


1.037530 


1.038064 


1.038597 


1.039128 


.06 


1.048267 


1.049105 


1.049924 


1.050734 


1.051537 


1.052327 


1.053108 


1.053883 


1 . 054653 


1.055417 


1.056178 


.07 


1.065511 


1.066641 


1.067750 


1.068844 


1.069926 


1.070995 


1.072054 


1.073103 


1.074146 


1.075181 


1.076211 


.08 


1.085296 


1.086759 


1.088198 


1.089613 


1.091015 


1.092401 


1.093775 


1.095137 


L. 096491 


1.097837 


1.099177 


.09 


1.107583 


1.109415 


1.111219 


1.112996 


1.114756 


1.116497 


1.118223 


1.119937 


1.121641 


1 . 1 23337 


1.125025 


.10 


1.132330 


1.134568 


1.136772 


1.138949 


1.141103 


1.143235 


1.145352 


1. 147455 


1.149547 


1.151630 


1.153705 


.11 


1.159496 


1.162177 


1.164820 


1.167429 


1.170012 


1.172573 


1.175115 


1.177643 


1.180161 


1.182668 


1.185170 


.12 


1.189048 


1.192205 


1.195318 


1 . 1 98394 


1.201441 


1.204465 


1.207470 


1.210461 


1.213439 


1.210110 


1.219376 


.13 


1.220947 


1.224613 


1.228230 


1.231807 


1.235353 


1.238874 


1.242376 


1.245865 


1.249343 


1 . 25281 6 


1 . 256285 


.14 


1.255166 


1 . 259371 


1.263523 


1.267633 


1.271711 


1.275764 


1.279800 


1.283822 


1.287838 


1.291851 


1.295863 


.15 


1.291674 


1.296449 


1.301168 


1.305844 


1.310487 


1.315106 


1.319710 


1.324305 


1 . 328896 


1.333488 


1 . 338084 


.16 


1.330448 


1.335824 


1.341140 


1.346414 


1.351656 


1.356876 


1.362084 


1.367288 


1 . 372493 


1.377705 


1.382928 


.17 


1.371471 


1.377477 


1.383422 


1.389325 


1.395199 


1 . 401055 


1.406905 


1.412756 


1.418616 


1.424490 


1.430382 


.18 


1.414728 


1.421393 


1.427999 


1.434563 


1.441104 


1.447633 


1.454162 


1 . 460700 


1.467255 


1.473834 


1.480442 


.19 


1.460212 


1.467565 


1.474861 


1.482122 


1.489364 


1.496602 


1.503850 


1.511116 


1.518411 


1 . 525741 


1.533111 


.20 


1.507918 


1.515991 


1.524010 


1.532000 


1.539980 


1 . 547966 


1.555973 


1.564010 


1.572090 


1.580217 


1.588400 


.21 


1.557849 


1.566672 


1.575448 


1.584202 


1 . 592958 


1 . 001732 


1.610540 


1.619395 


1.628306 


1.637282 


1 . 646330 


.22 


1.610014 


1.619619 


1.6)29185 


1.638741 


1.648312 


1.657917 


1.667572 


1.677290 


1 . 687084 


1 . 696961 


L. 706930 


.23 


1.664427 


1.674847 


1.685239 


1 . 695636 


1 . 706063 


1.716543 


1 . 727092 


1.737726 


1.748455 


1.759290 


1.770238 


.24 


1.721106 


1.732377 


1.743634 


1.754912 


1.766240 


1 . 777642 


1.789136 


1.800738 


1.812460 


1.824312 


1. 836)303 


.25 


1.780080 


1.792238 


1.804400 


1.816603 


1.828880 


1.841254 


1.853748 


1.866375 


1.879151 


1.892085 


1.905185 


.26 


1.841381 


1.854465 


1.867575 


1.880751 


1.894027 


1 . 907428 


1.920079 


1 . 934694 


1.948588 


1 . 962672 


1 . 976953 


.27 


1.905048 


1.919102 


1.933207 


1.947405 


1.961734 


1.976222 


1.990802 


2.005761 


2. 020844 


2.036152 


2.051602 


.28 


1.971129 


1.986199 


2.001348 


2.016625 


2.032066 


2. 047703 


2. 063560 


2. 079665 


2. 096003 


2.112614 


2.129495 


.29 


2. 039079 


2.055815 


2. 072063 


2. 088477 


2. 105095 


2. 121950 


2. 139068 


2.156467 


2. 174161 


2.192161 


2. 210474 


.30 


2.110759 


2.128016 


2. 145425 


2. 163041 


2.180906 


2.199055 


2.217513 


2. 236300 


2. 255430 


2.274913 


2. 294753 


.31 


2.184442 


2. 202881 


2.221516 


2. 240405 


2. 259595 


2.279119 


2. 299006 


2.319274 


2. 339937 


2.361003 


2. 382477 


.32 


2. 260809 


2. 28049(i 


2. 300430 


2. 320672 


2.341270 


2. 362262 


2. 383673 


2. 405523 


2. 427824 


2.450585 


2.473806 


.33 


2. 339949 


2. 360959 


2. 382273 


2. 403955 


2. 426056 


2.448613 


2. 471654 


2.495197 


2.519253 


2. 54382X 


2. 568923 


.34 


2. 421965 


2. 444380 


2. 467163 


2. 490383 


2. 514090 


2. 538324 


2.563111 


2. 588468 


2. 614407 


2. 640929 


2. 668033 


.35 


2. 506968 


2. 530881 


2. 555234 


2. 580098 


2. 605527 


2.631559 


2. 658220 


2. 685528 


2. 713488 


2.742103 


2. 771366 


.36 


2. 595085 


2. 620598 


2.646633 


2. 673262 


2. 700540 


2. 728506 


2. 757184 


2. 786589 


2.816726 


2. 847593 


2.879180 



Table 3. The correction factor for the second overtone of flexural vibration, T 3 , as a function of the diameter to length ratio, d/l, 

and Poisson's ratio, jx. 



The values of T 3 for d//<0.04 were obtained by a graphical interpolation as explained in the text. All other values were obtained from the computer solution. 


rf/z\u 


0.00 
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0.15 
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0.30 


0.35 


0.40 


0.45 


0.50 


0.00 


1.000000 


1. 000000 


1. 000000 


1. 000000 


1.000000 


1. 000000 


1. 000000 


I. 000(101) 


1. 000000 


1. 000000 


1. 000000 


.01 


1. 002600 


1.002000 


1.002710 


1. 002759 


1.002808 


1.002856 


1.002903 


1.002950 


1.002097 


1. 003043 


1.003089 


.02 


1.010421 


1. 010623 


1.010823 


1.011020 


1.011214 


1.011405 


1.011594 


1.011780 


1.011968 


1.012152 


1. 012335 


.03 


L. 023390 


1. 023845 


1. 024290 


1. 024729 


1.025164 


1. 025593 


1.026015 


1.026436 


1. 026855 


1. 027267 


1. 027676 


.04 


1. 041446 


1. 042248 


1.043034 


1. 043810 


1. 044578 


1. 045334 


1. 046082 


1. 046824 


1. 047563 


1. 048293 


1. 049019 


.05 


1. 064495 


1.005736 


1. 066954 


1. 068155 


1. 069343 


1. 070516 


1. 071678 


1. 072830 


1. 073974 


1.075110 


1. 076240 


.06 


1. 092436 


1. 094201 


1. 095938 


1. 097648 


1. 099340 


1. 101014 


1. 102673 


1. 104318 


1. 105953 


1. 107578 


1. 109196 


.07 


1. 125153 


1. 127523 


1. 129857 


1. 132160 


1. 134438 


1. 136693 


1. 138930 


1. 141151 


1. 143360 


1. 145557 


1. 147746 


.08 


1. 162526 


1. 165581 


1. 168592 


1. 171565 


1. 174505 


1. 177421 


1. 180314 


1. 183191 


1. 186053 


1. 188904 


1. 191747 


.09 


1. 204448 


1. 208262 


1. 212021 


1. 215737 


1.219416 


1. 223067 


1. 226696 


1. 230306 


1.233902 


1. 237488 


1. 241067 


.10 


1. 250814 


1. 255455 


1. 260034 


1. 264564 


1. 269056 


1. 273517 


1. 277955 


1. 282376 


1. 286786 


1.291188 


1. 295587 


.11 


1.301530 


1.307065 


1. 312532 


1.317947 


1. 323321 


1. 328666 


1.333000 


1. 339300 


1. 344603 


1. 349904 


I.35520S 


.12 


1.356516 


1. 363010 


1. 369431 


1. 375799 


1. 382127 


1. 388429 


1. 394714 


1. 400992 


1. 407270 


1. 413554 


1. 419850 


.13 


1. 415707 


1. 423222 


1. 430664 


1. 438052 


1. 445405 


1. 452737 


1. 460061 


1. 467387 


1. 474723 


1. 482077 


1. 489456 


.14 


1. 479051 


1. 487650 


1.496175 


1. 504653 


1.513102 


1. 521541 


1. 529982 


1. 538439 


1. 546921 


1. 555435 


1. 563990 


.15 


1. 546513 


1. 556257 


1.565933 


1. 575568 


1. 585187 


1. 594809 


1. 604449 


1. 614123 


1. 623839 


1. 633609 


1. 643439 


.16 


1. 618072 


1. 629024 


1.639915 


1. 650780 


1. 661643 


1. 672528 


1. 683453 


1.694432 


1. 705478 


1. 716601 


1. 727810 


.17 


1. 693723 


1. 705945 


1. 718120 


1. 730286 


1. 742473 


1. 754704 


1. 767001 


1. 779381 


1. 791856 


1. 804437 


1.817134 


.18 


1. 773475 


1. 787033 


1. 800562 


1.814105 


1. 827696 


1. 841361 


1. 855124 


1. 869002 


1.883011 


1. 897161 


1.911461 


.19 


1. 857353 


1. 872312 


1.887269 


1. 902268 


1. 917349 


1. 932540 


1.947868 


1. 963350 


1. 979003 


1. 994838 


2. 010865 


.20 


1. 945396 


1. 961827 


1. 978286 


1 . 994826 


2.011487 


2. 028302 


2. 045298 


2. 062497 


2. 079912 


2. 097557 


2. 115440 


.21 


2. 037655 


2. 055631 


2. 073675 


2. 091843 


2.110181 


2. 128725 


2. 147502 


2. 166535 


2. 185839 


2. 205426 


2. 225305 


.22 


2. 134200 


2. 153799 


2. 173513 


2. 193404 


2. 213523 


2. 233906 


2. 254584 


2. 275579 


2. 296906 


2. 318578 


2. 340601 


.23 


2.235112 


2. 256416 


2. 277894 


2.299610 


2.321620 


2. 343962 


2. 366670 


2. 389764 


2. 413260 


2. 437169 


2. 461496 


.24 


2. 340488 


2. 363589 


2. 386928 


2. 410580 


2. 434600 


2. 459032 


2. 483908 


2. 509250 


2. 535071 


2. 561381 


2. 588182 


.25 


2. 450441 


2. 475435 


2. 500745 


2. 526451 


2. 552613 


2. 579275 


2. 606470 


2. 634219 


2. 662535 


2. 691423 


2. 720883 


.26 


2. 565101 


2. 592092 


2. 619491 


2. 647382 


2. 675827 


2. 704873 


2. 734552 


2. 764882 


2. 795876 


2. 827534 


2. 859853 
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Table 4. A comparison of theoretical and experimental values 
of the correction factors Ti, T2, and T3 for the first 3 modes of 
flexural vibration of cylinders 





Theoretical 


Experi- 
mental 


Percent 
difference 


Specimen A, rf/f=0.27871_„ 
Specimen B, tf/Z=0. 45039. __ 


Tx 


1. 3787 
2. 0503 
2. 9814 

1. 9257 


1.3790 
2. 0498 
2. 9770 

1.9251 


-0.02 
+.02 
-.15 

+.03 



The specimens were cut from a single bar of soft steel, and had a value of 
Poisson's ratio of 0.2879. 

obtain and is probably in error by a larger amount 
than the other values. Also, it was necessary to 
extrapolate beyond the range of the tables to obtain 
this theoretical value. In spite of this, the indi- 
cated error is only 0.15 percent. 

We may conclude from this evidence that no 
significant disagreement with experiment can be de- 
tected by the use of presently available techniques. 
It is therefore believed that the estimate of error 
given herein is realistic, not only with regard to the 
errors in the numerical computations, but also with 
regard to the equations used. 



The author expresses his appreciation to Irene 
Stegun and Ruth Zucker for writing the IBM 704 
computer program used in these computations. 
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Infrared spectrum of hydrobromic acid, E. K. Plyler, J. Re- 
search NBS 64A, No. 5, 377 {I960). 

A precise measurement of the infrared fundamental bands of 
HBr 79 and HBr 81 has been made. The two band centers 
have been found to be 2558.94 and 2558.56 cm -1 , respectively. 
Rotational and vibrational constants have been calculated 
from the observed data. The constants are in good agree- 
ment with previous reported values. The centers of the two 
harmonic bands were used to calculate o> e x e and co e y e and they 
were found to be 45.58 and 0.072 cm -1 , respectively for 
HBr 79 and 45.56 and 0.072 cm" 1 for HBr 81 . 

The spectrum of singly ionized atomic iodine (I n), W. C. 

Martin and C. H. Corliss, J. Research NBS 64 A, No. 6, 443 
(1960). 

The I ii spectrum has been excited in electrodeless lamps 
and photographed from 655 to 11084 A. Wavelengths and 
estimated intensities are given for almost 2400 lines. A 
revision and extension of the earlier analyses of this spectrum 
has increased the number of known even levels from 43 to 
124, and the number of odd levels from 55 to 190. New 
gj-factors are given for 46 levels, and the previous designa- 
tions of 40 levels are changed. Improved measurements in 
the vacuum ultraviolet region given a correction of 7.4 cm -1 
to be subtracted from the values of all higher levels listed in 
Atomic Energy Levels, Vol. 3, (1958). The approximately 
1800 classified lines now include all of the strongest lines. 
The ^o of the ground configuration 5s 2 5p 4 has been found, 
and this configuration has been fitted to intermediate cou- 
pling theory. Magnetic dipole transitions between levels of 
the ground configuration, *F 2 - l D 2 (7282 A) and sp^So 
(4460 A), have been observed and their nature confirmed by 
the Zeeman effect. The line 5p 4 3 P 2 - 1 D 2 shows hyperfine 
structure which is in approximate agreement with a theo- 
retical calculation of the expected structure. New levels have 
been found for almost all higher configurations. All previ- 
ously known series have been extended and new ones found. 
From one of the new series, 5p 3 ( 4 S°)5-12(/ 5 G° the principal 
ionization energy for I u (154304 ± 1, cm -1 ) has been de- 
rived. The results of the analysis are compared with theo- 
retical expectations in a number of cases. 

The third spectrum of gold (Au m), L. Iglesias, /. Research 
NBS 64 A, No. 6, 481 {I960). 

The spark spectrum of gold has been photographed in a helium 
atmosphere from 500 to 6600 A. About 500 lines have been 
assigned to the third spectrum, Au in, and separated from 
those belonging to other stages of ionization, by observation 
of the polarity of the lines. Sixty two levels have been 
found: 17 even levels, arising from the 5d 9 and 5d 8 6s con- 
figurations; and 45 odd levels, belonging to the 5d 8 6p and 
5d 7 6s 6p configurations. All of the expected levels from the 
configurations 5d?, 5d 8 6s and 5d 8 6p have been identified 
except for the very high terms based on the 5d s (*S) core of 
Au iv. With these levels it was possible to classify 256 lines. 

Note on particle velocity in collisions between liquid drops 
and solids, O. G. Engel, J. Research NBS 64A, No. 6, 497 {I960) . 
Equations are developed for plane-wave particle velocity 
produced in solid-against-liquid collisions. An explicit 
expression for the dimensionless coefficient a that appears 
in these equations is deduced. 

Error analysis of a standard microwave phase shifter, G. E. 

Schafer and E. W. Beatty, /. Research NBS 64C, No. 4, 

261 {I960). 

A standard microwave phase shifter has been proposed 
which utilizes an adjustable short circuit attached to a 
tunable three-arm waveguide junction. Two classes of 



errors considered are those due to imperfect tuning of the 
waveguide junction and those due to inaccuracies in deter- 
mining the motion of the adjustable short circuit •and the 
width of the waveguide, termed tuning and dimensional 
errors, respectively. Graphs are presented for estimating 
the scattering coefficient magnitudes needed to estimate 
limits of tuning error from observations of amplitude changes 
during the tuning procedure. Graphs are presented for esti- 
mating limits of dimensional error for WR-90 waveguide in 
the recommended frequency range of 8.2 to 12.4 kMc/s from 
tolerances of the axial motion of the short circuit and broad 
dimensions of the waveguide. 

Electron resistivity studies on the Athabasca Glacier, Alberta, 
Canada, G. V. Keller and F. C. Frischknecht, J. Research 
NBS 64D, No. 5, 439 {I960). 

The use of electrical methods for measuring ice thickness and 
properties on the Athabasca Glacier, Alberta, Canada, has 
been studied by the U.S. Geological Survey. Two methods 
for measuring resistivity were tried: one, a conventional 
resistivity method in which current was introduced gal- 
vanically into the glacier through electrodes, and the other 
an electromagnetic method in which a wire loop laid on the 
ice was used to induce current flow. Results of the galvanic 
measurements showed large variations in the resistivity of 
the ice; in a surface layer several tens of feet thick the resistiv- 
ity is between 0.3 and 1.0 megohm-meters, and under this 
layer, the resistivity of the ice is more than 10 megohm- 
meters. The resistivity of the surface ice is determined by 
its water content rather than by molecular resonance loss. 
The ice had no effect on the mutual coupling measurements 
in the frequency range from 100 to 10,000 cycles per second. 
As a consequence the electromagnetic data could be inter- 
preted simply in terms of ice thickness and bedrock resistivity. 

Amplitude distribution for radio signals reflected by meteor 
trails, A. D. Wheelon, ,/. Research NBS 64D, No. 5, 440 {I960) . 
The probability distribution for the envelope of the received 
signal composed of reflections from many meteor trails is 
derived theoretically. Both the effects of numerous, small 
meteors and the residual reflections from infrequent, large 
meteors are treated simultaneously. For the particular ex- 
ample of exponential decay of initial spikes which are them- 
selves distributed as the inverse square of their amplitudes, 
we find that the probability that the composite residual 
signal amplitude exceeds a prescribed level r is given by 



P{R>r) = 



1 



r 1+ _i__T' 

L ^(^Q) 2 J 



This function behaves as a Rayleigh distribution for small 
amplitude margins r. For the larger, less likely amplitudes 
it agrees with the result predicted by elementary analysis of 
isolated meteor reflections. Possible refinements of these 
results are also discussed. 

Computation and measurement of the fading rate of moon- 
reflected UHF signals, S. J. Fricker, R. P. Ingalls, W. C. 
Mason, M. L. Stone, and D. W. Swift, J. Research NBS 64D, 
No. 5, 455 {I960). 

A method is described for predicting the fast fading rate of 
moon-reflected signals. It is based entirely upon considera- 
tions of the observer-moon positions and relative motions. 
Experimental results which are in good agreement with the 
computed fading rates have been obtained from a moon- 
reflection experiment at a frequency of 412 Mc/s. Some 
possible implications of this method of interpreting fading 
rates are given. 
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On the theory of wave propagation through a concentrically 
stratified troposphere with a smooth profile, H. Bremmer, 

,/. Research NBS 64D, No. 5, 467 (1960). 

The W.K.B. approximation for the solution of the height- 
gain differential equation for a curved stratified troposphere 
is discussed in detail. The approximation depends mainly on 
a variable ui(r) which can be interpreted as the height de- 
pendent contribution of the phase for a field solution obtained 
by separation of variables. An expansion of ui(r) with the 
aid of partial integrations leads to further approximations 
which facilitate the determination of the eigen values, and of 
the amplitudes of the modes connected with the propagation 
problem. The influence of the refractive-index profile, if as- 
sumed as smooth, then appears to be restricted to a depend- 
ence on the surface values of this index and of its gradient 
insofar as propagation over the ground is concerned. Further, 
all height effects of elevated antennas can be expressed in 
terms of the distance to the corresponding radio horizon. 
This results in simple relations between the fields connected 
with two different refractive-index profiles, provided both 
profile coincide near the earth's surface. 

Polarization and depression-angle dependence of radar ter- 
rain return, I. Katz and L. M. Spetner, /. Research NBS 64D, 
No. 5, 483 (1960). 

A study of recent experimental results on radar back scatter- 
ing from land and sea surfaces indicate: (a) the polarization 
dependence of the normalized radar cross section, o- , of 
ocean surfaces cannot be explained by the usual "interference 
phenomenon" and (b) there is a distinct difference in the form 
of the depression-angle dependence in that cr for "smooth" 
surfaces follows a negative exponential whereas o- for "rough" 
surfaces drop off as the sine of the depression angle. 

Methods of predicting the atmospheric bending of radio rays, 

B. R. Bean, G. D. Thaver, and B. A. Cahoon, /. Research 
NBS 64D, No. 5, 487 (1960). 

Three methods for predicting the bending of radio rays when 
the refractive index profile above the surface layer is unknown 
have been developed recently by the authors. These methods 
are: a statistical technique for refraction at high initial ele- 
vation angles, estimation of bending from an exponential 
model of atmospheric refractive index, and a modification of 
the exponential model to account for the heavily weighted 
effects of anomalous initial refractive index gradients at small 
initial elevation angles. Each model is dependent upon the 
value of the refractive index at ground level or, in the case of 
superrefraction, the additional knowledge of the refractive 
index gradient next to the earth's surface. Each method 
works best in a particular range of initial elevation angles or 
meteorological conditions. The height and angular ranges of 
application of each method are checked by comparison with 
values obtained from 77 diverse refractive index profiles 
representative of wide climatic variation. It is found that 
the use of the best of the three methods will always result in 
a prediction of the total atmospheric bending within 10 per- 
cent for initial elevation angles from zero to 10 mr and to 
within 4 percent for initial elevation angles greater than 17 
mr (~1 degree). 

Use of logarithmic frequency spacing in ionogram analysis, 

G. A. M. King, J. Research NBS 64D, No. 5, 501 (1960). 
The use of logarithmic frequency spacing in the reduction of 
ionograms to electron density profiles brings several advan- 
tages. Among them is the fact that, when computing factors 
for the analysis, one need not determine the group refractive 
index. Formulae involving only the phase refractive index 
are presented; for the ordinary component one exact and one 
approximate formula are given, while for the extraordinary 
component there is an approximate formula . . . valid over a 
wide range of geomagnetic latitudes. There is a brief dis- 
cussion of quasi-longitudinal approximations to the extraor- 
dinary phase refractive index. 

Guiding of whistlers in a homogeneous medium, R. L. Smith, 
J. Research NBS 64D, No. 5, 505 (1960). 

The velocity of energy flow of whistlers in a homogeneous 
medium is computed as a function of wave-normal angles. 
The maximum allowable cone of ray angles approaches 19 29' 



at very low frequencies, decreases with frequency to a mini- 
mum of 11° at a wave frequency of one-fifth the gyrofre- 
quency, then increases to 90° at the gyrofrequency. The 
velocity of energy flow departs markedly from the longi- 
tudinal value except at very low frequencies or very small 
wave-normal angles. 

Shielding of transient electromagnetic signals by a thin 
conducting sheet, N. R. Zitron, J. Research NBS 64D, 
No. 5, 563 (I960). 

The shielding effect of a thin, horizontal imperfectly con- 
ducting sheet against the transient field of a vertical magnetic 
dipole when excited by a ramp function is investigated. The 
results are calculated by taking Laplace transforms of the 
frequency spectrum functions for the steady-state problem. 
The response to the ramp function is calculated and the 
significance of the results in shielding against surges is 
discussed. 

Cylindrical antenna theory, R. H. Duncan and F. A. Hinchey, 
/. Research NBS 64D, No. 5, 569 (1960). 

A partial survey of cylindrical antenna theory pertaining to 
a tubular model with a narrow gap is presented. The survey 
includes discussion of the theories of Hallen, King and Middle- 
ton, Storm, and Zuhrt. A conceptual relation between theory 
and experiment is described. The latter part of the article is 
concerned with a new Fourier series solution of the Hallen 
equation. This solution is developed in such a way that the 
expansion coefficients are the unknowns of a system of linear 
equations. The elements of the coefficient matrix are given 
by a highly convergent series. Numerical results are given 
for half and full wavelength antennas with half length to 
radius ratios of 60 and 500. These results compare quite 
closely with those obtained from King-Middleton theory. 

Vibration-rotation structure in absorption bands for the 
calibration of spectrometers from 2 to 16 microns, E. K. 

Plvler, A. Danti, L. R. Blaine, and E. D. Tidwell, NBS 
Mono 16 (1960) 20 cents. 

Suitable bands of common gases have been tabulated and 
remeasured wherever necessary from 2 to 16 microns to 
obtain an accuracy of about 0.03 cm -1 throughout the region 
and to provide good calibrating points at frequent intervals. 
Some 600 rotation-vibration lines are illustrated in 20 spectro- 
grams and wavenumbers are listed in companion tables 
with considerable intercomparison with w T orthy data obtained 
in other laboratories. The absorption bands were remeasured 
or calibrated by using either a precisely graduated grating 
circle or standard atomic lines with the fringe system formed 
by a Fabry-Perot interferometer. Characteristic features 
of the individual bands are discussed briefly and references 
to other publications are given. The substances used for 
calibration include H 2 0, C0 2 , CO, HC1, HBr, NH 3 , C 2 H 2 , 
CH 4 , N 2 0, and polystyrene film. (Reprinted from the Journal 
of Research of the National Bureau of Standards — A. Physics 
and Chemistry, Vol. 64A , No. 1, January-February 1960.) 

Dynamic measurements of the magnetoelastic properties 
of ferrites, V. E. Bottom, NBS Tech. Note 49 (PB151408) 
(1960) $1.00. 

The relations between the mechanical and magnetic proper- 
ties of a ferrite are derived for the small signal or reversible 
condition using thermodynamic principles. The equation of 
motion of a ferrite ring driven in its fundamental mode is 
set up and solved leading to the equivalent electrical circuit 
of the magnetostriction resonator. From mechanical meas- 
urements of the density and dimensions of the ring and 
electrical and frequency measurements of the resonator, the 
elastic modulus, permeability, magnetostriction coefficients 
and loss factors in the ferrite can be determined. Apparatus 
is proposed for performing the above measurements. 

Magnetic drum directory and programming system for 
codesorting letter mail, P. C. Tosini, NBS Tech. Note 50 

(PB151409) < (1960) $1.75. 

This report is an analysis and extension of the Rabinow 
Engineering Company's proposal for a magnetic drum file- 
directory and special purpose computer to be used in a test 
coded-address ("codesorting") letter-mail sorting system 
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scheduled for late 1960 installation at the Washington, D.C. 
Post Office. The report is composed of 6 major sections that 
successively examine the major considerations involved. A 
short introduction is presented first. Then the general 
problem of coding the fields composing addresses is discussed. 
This area is analyzed in terms of facilitating both the human 
coding process and the computer coded address-to-bin- 
number translation. Addresses are classified into two types: 
'standard" and "non-standard". When coded, their fields 
must be unique and yet be sufficiently flexible to specify a 
variety of addresses. They must furthermore adequately 
specify some more common forms of misaddressed mail that 
is presently correctly sorted manually. A slightly amended 
form of Coding Plan E-l is then presented and analyzed in 
terms of the conditions specified for coded addresses in 
general. 

After the presentation of the coding area, the local mail sort 
translation program of the Rabinow proposal is discussed 
in detail. It is then extended to include the outgoing mail 
translation program and non-standard and special address 
forms. The information and instruction characters neces- 
sitated by the programs are listed in detail and the effects 
of the file-directory source document preparation determined. 
A coding keyboard is then suggested. 

The programming is then reexamined and it is determined 
that all addresses can be treated as special cases of the 
standard address forms at an increased cost in memory space 
and with a saving in the special inst ruction symbols required. 
Finally, alternative methods of file-directory information 
storage are explored and evaluated in terms of memory space 
requirements, table look-up and directory access times, 
and conceptual simplicity. 

Carrier frequency dependence of the basic transmission loss 
in tropospheric forward scatter propagation, K. A. Norton, 
NBS Tech. Note 53 {PB161554) {I960) $1.00. 
A further interpretation is given of certain Lincoln Laboratory 
data obtained in an experiment using scaled antennas as 
presented in a recent letter to the Proceedings of the I.R.E. 
from Bolgiano. This paper has four objectives: first, to 
clarify the significance of these data from the standpoint of 
the engineer developing long-range tropospheric scatter 
systems; second, to apply a further statistical analysis to these 
data; third, to consider their significance as regards the theory 
of radio propagation through a turbulent atmosphere; and 
fourth, to describe a suitable method for the measurement of 
the meteorological parameters entering the theory. Based 
on this analysis of the Lincoln Laboratory data, it is concluded 
that the carrier frequency dependence of the basic transmission 
loss cannot be variable from hour to hour. 

Determination of a general index of effort in sorting mail by 
conventional methods, S. Henig, NBS Tech. Note 54 
{PB161555) {I960) 50 cents. 

The conventional method of sorting letters in stages is 
described. It is shown, through reference to available 
statistics, that about 90% of letters originating within the 
areas served by several large post offices will be involved in 
no more than five stages to be completely sorted to a carrier, 
firm or section of boxes. A method is presented for the 
accounting of each stage's individual letter readings which are 
defined as the units of measurement for the index of sorting 
effort. It is found that the upper bound average number of 
readings per letter for three large post offices varies from 3.16 
to 2.98 and that these bounds are applicable to 90.2 to 92.4% 
of the local originating first class letters. It is anticipated 
that the values of such indexes for all large post offices will be 
consistently close to the range already determined. 

Emission spectra of N 2 , 2 , and NO molecules trapped in 
solid matrices, H. P. Broida and M. Peyron, J. Chan. Phys, 
32, No. 4, 1068 {I960). 

Molecular systems are observed in emission in solid products 
from a gas discharge trapped at liquid helium temperature. 
Previous tentative molecular assignments have been checked 
with the help of isotopic substitutions of oxygen and nitrogen. 
The Herzberg system {A bands) of oxygen (A 3 S M + — X 3 2 g ~) 
is analyzed and the molecular constants are derived for a 
molecule trapped in a nitrogen matrix. Another system {M 
bands) is attributed to the NO molecule ( 4 II — A 2 II). 



The components of power appearing in the harmonic analysis 
of a stationary process, M. M. Siddiqui, presented at Symp. on 
Statistical Methods in Radio Wave Propagation, University of 
California, Los Angeles, June 18-20, 1958. Reprinted from 
"Statistical Methods in Radio Wave Propagation,^ {Pergamon 
Press, New York, N.Y., 1960). 

Suppose t hat from a continuous record of a stationary process 
over a bimeinterval T'seconds, 2n f- 1 equally spaced readings 
are taken and the usual harmonic analysis is performed. 
Assuming the true frequencies present in the process to be 
. . . , l/(:rn, 1/(27^, 1/7 7 , 2/7\ . . . , cycles per second, 
it is shown that the power ascribed to the frequency j/T c/s 
consists of three components: (1) the true power in the 
frequency j/T c/s; (2) the powers in the frequencies {2kn±j)/ 
T c/s, k=l,2, . . . ; (3) a part of the powers in the frequencies 
l/{mT) c/s, m=2, 3, . . . . 

Introduction — Can you measure it? L. V. Judson, Am. Soc. 
Toolcfc Mfgr. Engrs. Tech. Paper 239, Hook 1, 60, 1 (I960). 
A paper wmich is the introductory paper in a Metrology 
Symposium. It indicates the need for precise and exact 
measurements and it treats two subsidiary questions: (1) 
Why measure? and (2) What lias been the history of measure- 
ment? 

Correction for systematic wavelength shift in atomic beam 
devices, 11. L. Barger and K. G. Kessler, J. Opt. Soc. Am. 50, 
No. 4, 352 {I960). 

Light emitted or absorbed by atoms in an atomic beam is 
shifted in frequency relative to the frequency for a stationary 
atom if the light ray is not normal to the trajectory of the 
atom. The magnitude of this shift is calculated and several 
devices wdiich correct for this effect are described. For 
photoelectric recording of interference fringes, complete 
compensation can be achieved. In the case of photographic 
recording, it is shown that the residual error is negligible. 

The error in prediction of F2 maximum usable frequencies 
by world maps based on sunspot number, E. L. Crow and 
I). H. Zacharisen, Section II, Radio Propagation Phenomenolo- 
</!/, hook, Statistical methods in radio wave propagation, p. 248 
(Pergamon Press, New York, N. Y., 1960). 
The feasibility of preparing semi-permanent world maps for 
predicting F2 maximum usable frequencies based on sunspot 
number is studied. The components of variance of I he predic- 
tion are estimated; they arise from the following four major 
sources of error: determining regressions of P ,i characteristics 
on sunspot number with limited years of data, the scatter of 
points about these regressions, drawing world contour maps 
with limited numbers of stations, predicting the future sunspot 
number. All four sources are found to contribute about 
equally to a total standard deviation of monthly median 
maximum usable frequency predicted five months in advance 
of about 1 and 3 Mc for and 4000 km respectively. Hence 
the proposed maps are concluded to be feasible. 

Bremsstrahlung linear polarization, J. W. Motz and R. C. 
Placious, II Nuovo Cimento 15, Series X, 571 {I960). 
This work presents a general quantitative descrption of 
the bremsstrahlung linear polarization on the basis of original 
experimental data and available theoretical calculations. 
The results give the dependence of the polarization on a) the 
initial electron kinetic energy, r T in a range from 10~ 2 to 
10 3 MeV, b) the photon energy in a range from 0.1 T to 7'o, 
c) the photon emission angle in a range from to 180 degrees, 
and d) the atomic number of the target in a range from 4 to 
79. The experimental data were obtained for a range of 
electron energies from 0.05 to 1.0 MeV with beryllium, 
aluminum, and gold targets. Theoretical estimates of the 
polarization were calculated from the Sommerfeld-Kirk- 
patrick-Wiedmann results for the non-relativistic energy 
region, the Olsen-Maximon results for the extreme relativist ic 
energy region, and the Gluckstern-Hull (Born approximation) 
results for the intermediate energy region. Final results are 
expressed in terms of the peak polarization and the corre- 
sponding peak angle as a function of the electron and photon 
energies, and best estimates of the polarization are given on 
the basis of the combined experimental and theroretical data. 
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Geomagnetic disturbances and velocity of slow-drift solar 
radio bursts, M. G. Wood and C. S. Warwick, Nature 184, 
No. 4697, 1471 (1959). 

Frequency drift rates were measured on records of the Fort 
Davis dynamic spectrum analyzer for 48 Type II (slow-drift) 
solar radio noise bursts. The change with frequency of the 
frequency drift rate is found to differ for two groups of bursts : 
those followed by geomagnetic disturbance and those followed 
by no disturbance. This change in drift rate is interpreted as 
acceleration of the radio source as it moves outward through 
the solar corona. Bursts followed by geomagnetic disturb- 
ance determine a rate that corresponds to greater positive 
acceleration of the source. 



The outlook for machine translation, F. L. Alt, Proc. Western 
Joint Computer Conf., Volume 17, p. 203, San Francisco, 
Calif., (May 1960). 

The current status of the major research projects concerned 
with the use of computing machines for the translation of 
natural languages (e.g., Russian to English) is surveyed. 
None of these projects has as yet reached the stage of fully 
automatic translation of satisfactory quality. The difficulties 
encountered fall into two main classes: semantic and syn- 
tactic. Several methods of approach have been used to 
overcome them; these are surveyed, and the approach of the 
National Bureau of Standards is described in more detail. 
It appears likely that within a few years automatic translation 
of fair quality may be possible, at a cost substantially below 
that of human translation. The major cost elements — initial 
development of machine programs, cost of machines and of 
operation — are analyzed in their dependence on the method 
of translation and on the type of equipment used. 

Statistical evaluation of interlaboratory cement tests, J. R. 

Crandall and R. L. Blaine, Am. Soc. Testing Materials Proc. 
59, 1129 (1959). 

A cooperative series of physical and chemical tests was made 
by 103 laboratories on twelve samples of cement during the 
period of one year. Results of tests by all laboratories on 
each property were plotted in scatter-diagrams according to 
the Youden method to enable the participating laboratories 
to quickly evaluate their results. The results were treated 
statistically. Both statistical and graphic methods were 
employed to indicate the precision of the test methods and the 
extent of laboratory bias. The lack of precision with some 
of the tests indicated a need for improved test methods. 
The work of each of the various laboratories was evaluated 
by a rating system. A large number of laboratories obtained 
results that were in good agreement but a few laboratories 
showed poor agreement with the larger group. Many of the 
discrepant laboratories did not improve during the one-year 
period the test program was active. 

Discrepant results obtained by a few laboratories inflated 
appreciably the standard deviations among laboratories for 
the various tests. 

The need for a continuing reference sample program was 
indicated. 



Effect of notch geometry on tensile properties of annealed 
titanium at 100°, 25°, -78° and -196° C, G. W. Geil and 
N. L. Carwile, Am. Soc. Testing Materials Proc. 59, 985 (1959). 
The combined effects of low temperatures and geometry of 
a circumferential notch on the tensile behavior of initially 
annealed commercially-pure titanium was investigated. The 
notch geometry was varied by changes in its depth and root 
radius; the notch depth, in terms of the cross-sectional area 
removed in machining the notch, ranged from 5 to 87 percent 
and the root radius ranged from 0.005 to 2.0 in. True stress- 
true strain data were obtained to the initiation of the fracture 
of specimens at 100°, 25°, -78° and -196° C and repre- 
sentative true stress-true strain curves are presented. The 
relationships of selected strength and ductility indices to the 
notch geometry and to the accompanying triaxialities and 
stress concentration factors are presented graphically and 
discussed in some detail in the text. 



The wavefront reversing interferometer, J. B. Saunders, 
Colloquium on Optics and Metrology, Sponsored by the Intern, 
Commission of Optics, Brussels, Belgium, May 8—9, 1958. 
Publ. in Optics in Metrology, p. 227 (Pergamon Press, New 
York, N. Y., 1960). 

The wavefront reversing interferometer is designed to test 
the shape of wavefronts. If a wavefront is characteristic of 
the shape of an optical element, or combination of elements, 
the evaluation of the wavefront is indirectly an evaluation of 
the element. It serves to measure the shapes of reflecting 
surfaces, the aberrations of lenses, and the variations of 
optical paths in a beam, such as inhomogeneity of the medium 
through which the light traveled. A unique method of 
analyzing the data is given. An instrument suitable for 
testing lenses and mirrors, by projecting the interference 
fringes onto a screen so that they may be observed by an 
audience, is described. 

Spectra emitted from rare gas-oxygen solids during electron 
bombardment, L. J. Schoen and H. P. Broida, J. Chem. Phys. 
32, No. 4, 1184 (I960). 

Light emission from rare gas solids containing small amounts 
of oxygen has been excited by an electron bombardment 
technique. The Herzberg bands of oxygen and a line group 
believed to be the 1 S-^ 1 D transition of the neutral oxygen 
atom form the strongest features of the observed spectra. 
In a neon matrix, the Second Negative system of 02 + has 
been identified. The effect of the various solids on the 
vibrational structure of the Herzberg bands is discussed. 

Accurate microwave wavemeters with convenient calibration 
tables, H. E Bussey and A. J. Estin, Rev. Set. Instr. 31, 
410 (1960). 

Accurate and convenient microwave cavity wavemeters are 
described that are suitable for many precise physical measure- 
ments. High Q's were attained by refined construction 
techniques. At 9000 Mc a precision of 0.02 Mc was obtained. 
The absolute accuracy also may be very high after strains in 
the metal have stabilized. A calibration table containing 
10 4 entries, easily formed by means of a high speed computer, 
makes the wavemeters convenient to use. The curve fitting 
method, accurate to one in 10 6 , is described. 

The Stokes flow problem for a class of axially symmetric 
bodies, L. E. Payne and W. H. Pell, /. Fluid Mech. 7, PL 4, 

529 (1960). 

The Stokes flow problem is concerned with fluid motion about 
an obstacle when the motion is such that inertial effects can 
be neglected. This problem is considered here for the case 
in which the obstacle (or configuration of obstacles) has an 
axis of symmetry, and the flow at distant points is uniform 
and parallel to this axis. The differential equation for the 
stream function \p then assumes the form L 2 -i\p=0, where 
L-i is the operator which occurs in axially symmetric flows 
of the incompressible ideal fluid. This is a particular case 
of the fundamental operator of A. Weinstein's generalized 
axially symmetric potential theory. Using the results of this 
theory and theorems regarding representations of the solu- 
tions of repeated operator equations, the authors (1) give 
a general expression for the drag of an axially symmetric 
configuration in Stokes flow, and (2) indicate a procedure 
for the determination of the stream function. The stream 
function is found for the particular case of the lens-shaped 
body. 

Explicit calculation of the drag is difficult for the general 
lens, without recourse to numerical procedures, but is rela- 
tively easy in the case of the hemispherical cup. As ex- 
amples illustrative of their procedures, the authors briefly 
consider three Stokes flow problems whose solutions have 
been given previously. 

Analysis of two-factor classifications with respect to life tests, 

M. Zelen, Book, Contributions to Probability and Statistics, 
Paper 42, p. 508 (Stanford Univ. Press, Stanford, Calif., 1960). 
Recently, much attention has been given to a study of the 
exponential distribution with particular emphasis on appli- 
cations in life testing, e.g., Epstein and Sobel, Gupta, and 
Sobel. These papers have mainly considered techniques and 
procedures associated with observations from a single ex- 
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ponential distribution. This paper outlines procedures to be 
used for analyzing the results of a two-way classification with 
respect to the scale factor of the exponential distribution. 
The test procedures are based on likelihood ratio bests and 
have their analogs in the usual analysis-of-variance tests for 
main effects and interactions. Approximations are given for 
tests based on small samples. These same procedures carry 
oyer directly to a two-way analysis of variance, if the quan- 
tities to be analyzed are variance estimates following chi- 
square distributions. 

Microwave spectrum of trans-crotononitrile, V. W. Laurie, 
/. Chem. Phys. 32, No. 5, 1588 {I960). 

The microwave spectrum of ethyl cyanide has been studied 
in the region from 17-36 kMc. Both parallel and perpen- 
dicular transitions have been assigned. Rotational constants 
(Mc) for the ground vibrational state are a Q = 27 663.30, b Q = 
4714.14, c = 4235.14, £>., = 0.0035, D JK = -0.0496. From 
Stark effects the dipole moment is calculated to be yua = 3.78 D, 
fib=1.3S Z),;u=4.02 D. Hyperfine splittings lead to a value 
of —3.3 Mc for the quadrupole coupling of the N 14 nucleus 
along the a axis. Rotational transitions of molecules in 
several excited vibrational states have also been measured. 
A Coriolis interaction between the torsional mode and the 
in-plane CCN bend has been found and is discussed. Split- 
tings of transitions of molecules in the first excited torsional 
state show that the barrier hindering internal rotation of the 
methyl group is 3.05 kcal/mole. 

Disk distribution of flares associated with certain radio 
bursts, C. Warwick, Proc. Symp. Physical Processes in the 
Sun-Earth Environment, July 20-21, 1959 {Ottawa, Defence 
Research Telecommunications Establishment) , DRTE Pub No 
1025, p. 249 {I960). 

Flares associated in time with U-bursts and with slowdrift 
bursts have an unusual distribution over the solar disk that 
may be caused by a directive effect in the radio emission. 
Such an effect, making the observation of these radio bursts 
more probable from a source near the limb of the sun than 
from a source near the center, could arise in coupling processes 
either between plasma waves and electromagnetic radiation, 
or between the two modes of electromagnetic radiation. 

List of Titles 

Journal of Research, Section 64A, No. 5, September-October 
1960. 70 cents. 

Infrared spectrum of hydrobromic acid. E. K. Plyler. (See 

above abstracts.) 

Determination of the value of the faraday with a silver- 
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Systems silver iodide-sodium iodide and silver iodide-potas- 
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sorption spectra of some acetylated aldopyranosides. 
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Tolerances for layer thicknesses in dielectric multilayer 
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Heat of formation of decaborane. W. H. Johnson, M. V. 
Kilday, and E. J. Prosen. 

Ultra low-conductivity water by electrophoretic ion exclu- 
sion. W. Haller and H. C. Duecker. 

Spectrophotometric determination of the ionization constant 
of dimethylpicric acid (2,4,6-trinitro-3,5-xvlenol) in water 

^ at 25° C. M. M. Davis, M. Paabo, and R. A. Robinson. 

Spectrophotometric determination of the ionization constant 
of 2,4,6-trinitro-m-cresol in water at 25° C. M. M. Davis 
and M. Paabo. 

Method for the separation of titanium, zirconium, iron, and 
aluminum from one another and for their subsequent 
determination. T. J. Murphy, W. S. Clabaugh, and R. 
Gilchrist. 
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1960. 75 cents. 

Error analysis of a standard microwave phase shifter. G. E. 
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Journal of Research, Section 64D, No. 6, November-Decem- 
ber 1960. 75 cents. 

URSI National Committee Report: 

Commission 1. Radio Measurement Methods and Stand- 
ards 
Commission 2. Tropospheric Radio Propagation 
Commission 3. Ionospheric Radio Propagation 
Commission 4. Radio Noise of Terrestrial Origin 
Commission 5. Radio Astronomy 
Commission 6. Radio Waves and Circuits 
Commission 7. Radio Electronics 



Temperature-induced stresses in solids of elementarv shape, 
L. H. Adams and R. M. Waxier, NBS Mono. 2 (1960) 

2 § o f* ii i - ^ 

Table of wavenumbers, 2000 A to 7000 A, C. D. Coleman, 
W. R. Bozman, and W. F. Meggers, NBS Mono. 3, Vol. 

I (1960) $6.00. 

Table of wavenumbers, 7000 A to 1000 /*, C. D. Coleman, 
W. R. Bozman, and W. F. Meggers, NBS Mono. 3, Vol. 

II (1960) $6.00. 

Mercury barometers and manometers, W. G. Brombacher, 
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